Large tuneable Rashba spin splitting of a two-dimensional electron gas in Bi2Se3 
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We report a Rashba spin splitting of a two-dimensional electron gas in the topological insulator 
Bi2Se3 from angle-resolved photoemission spectroscopy. We further demonstrate its electrostatic 
control, and show that spin splittings can be achieved which are at least an order-of-magnitude 
larger than in other semiconductors. Together these results show promise for the miniaturization of 
spintronic devices to the nanoscale and their operation at room temperature. 
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Spintronics promises to revolutionize electronics and 
computing by making explicit use of the electron's spin in 
addition to its charge [T] . A key-requirement is spin ma- 
nipulation via an electric, rather than a magnetic, field. 
However, even the prototypical method [5] for demon- 
strating such control, via an electrostatically-tuneable 
Rashba [3] spin splitting of a two-dimensional electron 
gas (2DEG), is difficult to implement: despite decades 
of research, crucial spintronic components such as the 
spin-field and spin-Hall effect transistors have only been 
realized very recently in a laboratory setting [H |3] , and 
are currently restricted to operation at cryogenic tem- 
peratures. A breakthrough is hampered by the intrinsic 
properties of available materials. In particular, the mod- 
est spin-orbit interaction of most semiconductors leads to 
small Rashba splittings, necessitating low temperatures 
for device operation and long channel lengths with ultra- 
high purity material to avoid spin-flip scattering events. 

States with larger Rashba splittings are known to ex- 
ist at several metal surfaces |^{S] and in ultrathin metal 
films [HI HO]- Unlike a semiconductor 2DEG, however, 
their spin splitting cannot be influenced by external elec- 
tric fields, rendering such systems unsuitable for real de- 
vices. Here, we obtain a Rashba splitting of a semicon- 
ductor 2DEG at the surface of the topological insulator 
Bi2Se3. Through electrostatic control, we achieve spin- 
splittings comparable to those of metal surface states and 
surface alloys, and orders of magnitude larger than in 
other semiconductors. This system should therefore pro- 
vide an ideal platform to establish spintronic control in 
practical devices. 



In-situ cleaved single-crystal samples of Bi2Se3 and Ca- 
doped Bi2Se3 were measured using conventional, time- 
of- flight, and spin-resolved angle-resolved photoemission 
(ARPES). The time-of-flight measurements were per- 
formed using a VG Scienta ArTOF 10k spectrometer [TT] 
at BESSY II with the storage ring operating in single 
bunch mode at a repetition rate of 1.25 MHz. Gon- 
ventional ARPES measurements were performed using 
a SPEGS Phoibos 150 hemispherical analyser at the 
ASTRID synchrotron [T^. The spin- resolved measure- 
ments were performed using the GOPHEE setup at the 
Swiss Light Source [13]. The energy and angular resolu- 
tion was better than 5-15 meV (80 meV) and 0.13-0.2° 
(1.5°) for the spin-integrated (resolved) measurements, 
respectively. 

The freshly-cleaved sample (Fig. [T|a)) shows the fa- 
miliar linearly-dispersing topological surface state, whose 
apex defines the Dirac point, as well as the bulk valence 
bands at higher binding energy. Close to the Fermi level, 
occupied bulk conduction band states are also visible, 
due to the usual degenerate electron doping of this com- 
pound. With time, the Dirac point shifts to higher bind- 
ing energies (Fig. [ijb)), indicating a downwards shift of 
the electronic bands close to the surface, associated with 
an increased surface electron doping [TH US]. The re- 
sulting quantum well causes a 2DEG to form at the sur- 
face [14] . with two-dimensional subbands giving rise to 
the well-defined Fermi surface in Fig. [ijb). 

With further surface doping (see Supplementary Movie 
1 [16]), the band bottom of the 2DEG state starts to shift 
away from the F-point. This is indicative of a Rashba 
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FIG. 1: (a-c) Electronic structure of Bi2Se3, measured 
using time-of-flight ARPES at increasing times (marked in 
Fig.|4];a)) after cleaving the sample {hu = 19.2 eV, T = 10 K). 

(d) Spin-ARPES total intensity MDC at the Fermi level 
{hu = 19.5 eV, T = 60 K) along the f-M direction, and 

(e) the measured (points) and fitted (lines) radial {Px), tan- 
gential (Py), and out-of-plane (P^) components of the spin 
polarization along this cut. (f) Schematic representation of 
the in-plane component of the Fermi surface spin texture. 



spin splitting: in the presence of a potential gradient 
(provided here by the band bending) , spin-orbit coupling 
lifts the spin-degeneracy of the 2DEG, separating the 
states into two bands. In the simplest approximation, 
their dispersion is given by 

E^m)=Eo + —^±a^, (1) 
2m* 

where m* is the effective mass and a is the Rashba cou- 
pling parameter, dependent on both the gradient of the 
potential and the spin-orbit coupling strength. The spin- 
split bands cross at fcy = 0, causing the band bottom to 
shift away from the zone centre as observed here. As the 
band bending, and consequently the near-surface poten- 
tial gradient, continues to increase, a well defined spin 
splitting of the states is established (Fig. [TJ;), much like 
for the model example of the Au(lll) surface state [5]. 
We stress, however, that the states we observe here are 
not surface states |17j , but rather quantum-confined con- 
duction band states, that is, a 2DEG similar to those 
found at a number of semiconductor surfaces and in- 
terfaces [IHl m]. However, the size of the splitting is 
much larger. Even for semiconductors generally consid- 



ered to have a strong spin-orbit coupling, such as InAs 
and InSb, no Rashba splitting of a 2DEG could be re- 
solved by previous ARPES studies [20l[2T]. Indeed, the 
momentum (energy) splitting at the Fermi level of up 
to Akp = 0.08 A-i {AEr = 180 meV) are approxi- 
mately one to two orders of magnitude larger than in 
semiconductor 2DEGs which have previously been con- 
sidered good candidates for spintronic devices p2H24] . 
and a factor of 2-4 larger than for the Au(lll) surface 
state [3 [55]. 

Apart from expected matrix-element effects |14j . nei- 
ther the topological state nor the spin-split states vary 
when measured with different photon energies (Supple- 
mentary Movie 2 [E]). This lack of kz dispersion con- 
firms the two-dimensional nature of the electron gas that 
gives rise to the Rashba splitting. Further, we have con- 
firmed their spin-split nature by spin-resolved ARPES. 
The spin polarization of the bands crossing the Fermi 
level along the T-M direction (Fig. [ijd)) were extracted 
(Fig. [TJe)) using the procedure described in Ref. [26] . 
The total intensity momentum distribution curve (MDC) 
was fit with 7 components, accounting (from high |fc|||) 
for the topological state (gray lines) and the spin-split 
2DEG state (black lines), with a single central compo- 
nent to account for the residual bulk conduction band 
and unresolved shallow 2DEG states. The analysis shows 
that both the topological and the spin-split 2DEG states 
are 100% spin polarized and the spin texture of the lat- 
ter is that expected for a simple Rashba model. The 
spins lie predominantly in-plane, and averaging over the 
Kramers pair are situated at an angle of ~— 95° (~86°) 
relative to the state's momentum. The deviations from 
perfect spin-helicity most likely arise from a slight sample 
misalignment, although we cannot rule out spin-canting 
due to hexagonal warping of the Fermi surface. How- 
ever, a detailed analysis of such features is beyond our 
current experiment. The spin of the topological state is 
aligned approximately anti-parallel to that of the outer 
Rashba-split branch of the 2DEG state, and possesses fi- 
nite out-of-plane spin polarization, possibly [27] due to 
hexagonal warping of this state once the 2DEG starts to 
develop [T3] . The in- plane spin texture around the Fermi 
surface extracted from these measurements is summa- 
rized schematically in Fig. [ijf). 

Tuning matrix elements via the photon energy, we 
can resolve additional shallow states which are also non- 
dispersive in kz (Supplementary Movie 2 |16j). From 
a detailed measurement (Fig. [2]^ a)), a clear second state 
can be discerned above the spin-split band bottom, while 
a very weak and diffuse third state can just be distin- 
guished close to the bulk conduction band minimum. In 
agreement with earlier studies [Ml [201 HH] > we attribute 
these to a ladder of quantum- well subband states induced 
by the surface band bending. To support this assignment, 
we have performed Poisson-Schrodinger calculations, fol- 
lowing the method described in Refs. [UJ [55]. Taking 
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FIG. 2: (a) Detailed measurement {hu = 16 eV, T = 60 K) 
of multiple 2DEG states (upper features) and the topologi- 
cal surface state (lower feature). Model Poisson-Schrodinger 
calculations (red solid lines) reproduce the ladder of 2DEG 
states within the surface quantum well. Inset: EDC at T 
with 2DEG states marked by arrows, (b) Near-surface band 
bending (blue) with the calculated energy levels and modulus- 
squared wave functions of the quantum-confined states (red). 
A Rashba splitting of a = 0.36 eVA (estimated from the kp- 
splitting) applied to the lowest calculated subband dispersion 
gives the green lines in (a). 



the total band bending magnitude as the shift of the 
Dirac point from its position in the freshly cleaved sam- 
ple, and taking the same materials parameters detailed 
in Ref. |14j . we find quite good agreement vi^ith the ex- 
perimental data. This agreement is further improved by 
a small reduction of the static dielectric constant to 70, 
yielding the subband structure shown in Fig. [2] strongly 
supporting the multiple quantum- well-state origin of the 
observed bands. The change in dielectric constant com- 
pared to our earlier work |14j may result due to the in- 
creased free-carrier density within the 2DEG established 
here, a factor which is known to cause a large reduction 
of dielectric constant in other semiconductors |30| . 

Only the lowest of these measured subbands exhibits 
a visible spin splitting in the ARPES. In a simple model 
developed to describe Rashba splitting of III-V semicon- 
ductor quantum wells [31J, the Rashba coupling param- 
eter 
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e'{z) e'(z) + A 



VX^) ) , (2) 



where £'{z) ~ e + V{z)+Eg, e is the subband energy rela- 
tive to the bulk conduction band, V{z) is the band bend- 
ing potential. Eg the band gap, and A is a measure of 
the spin-orbit coupling within a Kane k • p approach [32] ■ 
Consequently, for the deepest subband, whose wave func- 
tion is localized close to the surface in the region of 
steepest potential (Fig. [2]jb)), a much stronger Rashba 
splitting would be expected than for the shallower states, 
whose wave functions are much more extended in regions 
of shallower potential. Indeed, taking a representative 
value of A = 1 eV, this model predicts Rashba parame- 
ters of 0.1, 0.035, and 0.014 eVA for the three subbands. 
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FIG. 3: Rashba-split 2DEG in Ca-doped (non-degenerate) 
Bi2Se3 induced by (a) 18 L of CO-adsorption at 60 K and (b) 
<C 1 ML of Rb-deposition at 170 K. The measurement temper- 
ature in (a) and (b) is 60 K and 380 K, that is, well above room 
temperature, respectively, (c-e) Reversible nature of effective 
gating of undoped Bi2Se3, with the spin splitting (along F- 
K) changing from (c) Afc^ = 0.034 A"^ [Vg ~ 0.28 eV) to 
(d) AfcF = 0.011 A~^ (Vg ~0.11 eV) following annealing the 
sample at 250 K to remove some surface adsorbates to (e) 
A^F = 0.025 A~^ as the effective gate potential is increased 
again to Vg ~0.26 eV. 



This underestimates the observed splitting of the lower 
subband by a factor of ~ 4, which is not unexpected for 
a simple model intended to describe standard semicon- 
ductor materials with quite different bulk band structure 
and much weaker spin-orbit interactions. However, the 
qualitative agreement with the trend of spin splittings 
between the different subbands strongly supports that 
the Rashba splitting here is dominated by the gradient 
of the back-side of the potential well. 

Most likely, the 2DEG formation here is mediated by 
the adsorption of residual gas molecules such as CO from 
the vacuum onto the surface. Indeed, deliberate dos- 
ing with CO leads to much more rapid formation of the 
2DEG [33] (Fig. [Sja)), while its formation can also be 
accelerated via dehberate adsorption of H2O [34], an- 
other common residual gas in ultra-high vacuum. The 
spin-split 2DEG can also be induced by deposition of 
small quantities (^ 1 ML) of alkali metals on the sur- 
face (Fig. [3]^b)), as well as by deposition of other donor 
species such as Fe [17] or Gu [35]. Irrespective of the ex- 
act microscopic origin of the surface donors, their pres- 
ence is analogous to application of an external gate field. 
Indeed, via removal of surface adsorbates by mild an- 
nealing (250 K), we reduce this effective gate field, show- 
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FIG. 4: (a) Evolution of the Fermi wave vector along f — M 
and T — K with time after cleaving the sample, (b) Extracted 
kp positions and band-bending which induces the 2DEG. (c) 
Sheet density and (d) fcy-splitting at the Fermi level of the 
measured (points) and calculated (lines) Rashba-split 2DEG 
as a function of the effective gate potential, Vg. 



ing that the spin spUtting can correspondingly be dimin- 
ished. This demonstrates the reversible and controllable 
nature of the Rashba splitting, driven by a change in the 
band bending (Fig. [sjjc-e)). We note that this is funda- 
mentally different from the adsorbate-induced changes to 
Rashba splittings achieved in conventional surface states, 
as here it reflects the electrostatically-tuneable nature of 
a semiconductor 2DEG. In a practical device, of course, 
this would be achieved by electrically-biasing a gate elec- 
trode. The injected charge necessary to induce the strong 
Rashba splittings observed here may be achievable us- 
ing conventional dielectric media, where substantial gate- 
tuned changes in Hall concentration have already been 
demonstrated for Bi2Se3 thin films [551. Alternatively, 
the sheet densities required are well within the range 
that have recently been realized in electric-double-layer 
FETs [37], which is directly analogous to the surface 
alkali- metal doping shown in Fig. . 

Using surface adsorbates to achieve this field-effect 
control uniquely permits the simultaneous spectro- 
scopic investigation of the 2DEG. Utilizing time-of- flight 



ARPES [TT], we track the continuous evolution of the 
full {kx,ky, E)-dependein.t electronic structure as a func- 
tion of effective gate voltage (Fig. |4] and Supplementary 
Movie 1 [H]). With time, the effective gate potential, 
given by the shift of the Dirac point, monotonically in- 
creases in magnitude up to ^ 0.35 eV. Concurrent with 
this, the Fermi surface of the 2DEG grows and smoothly 
evolves into two spin-split sheets above an effective gate 
voltage of ^0.15 eV. Crucially, through a combination of 
different sample temperatures, background pressure, and 
differing times after the cleave, the effective gate poten- 
tial reaches a factor of three higher than in our previous 
study [14], a necessary condition for realization of the 
Rashba split states we report here. From these measure- 
ments, the sheet density of the 2DEG (Fig.[4]Jc), obtained 
from its Luttinger area) and the momentum splitting at 
the Fermi level (Fig. |4|^d)) can be extracted as a function 
of the effective gate voltage. The measured sheet density 
and trend in the spin splittings are well reproduced by 
our model calculations, again supporting the importance 
of the potential gradient of the backside of the quantum 
well in determining the Rashba splitting. 

By tuning this gradient by an effective gate voltage 
change of only ~ 200 mV, the spin splitting can be varied 
from zero up to ~0.08 (Fig.|4|d), corresponding to a 
Rashba parameter of ^ 1.3 eVA. This is orders of magni- 
tude larger than the control of spin-splitting achieved in 
conventional semiconductor 2DEGs [23], making Bi2 Sea- 
based 2DEGs ideal for applications where electrical ma- 
nipulation of the spin precession is required. For exam- 
ple, the precessional phase shift for injected electrons in 
the spin-FET is AkpL, where L is the channel length [2]. 
Consequently, for the 2DEG created here, the required 
channel length for maximal source-drain current modu- 
lation could be as small as ~ 4 nm, providing a route 
towards nanoscale spintronic devices. Furthermore, the 
large Rashba effects ensure that a well defined spin- 
splitting is maintained up to well above room temper- 
ature, as shown in Fig. |3][b). This is a necessary re- 
quirement for practical implementation of spintronics, 
but something which has yet to be demonstrated for de- 
vices based on more conventional semiconductors [H [S] . 
The recent achievements of electrostatic gating [SH] [SS] 
of Bi2Se3 as well as its epitaxial growth on Si |39j and 
GaAs [40 indicate good prospects to incorporate the 
spintronic functionality suggested here into practical de- 
vices, and to integrate these with conventional semicon- 
ductor electronics. 
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